Zero-Temperature Properties of the Quantum Dimer Model on the Triangular Lattice 



Arnaud Ralko, 1 Michel Ferrero, 2 Federico Becca, 2 Dmitri Ivanov, 1 and Frederic Mila 1 
1 Institut de Theorie des Phenomenes Physiques, Ecole Polytechnique 
Federate de Lausanne (EPFL), CH-1015 Lausanne, Switzerland 
2 INFM-Democritos, National Simulation Centre and International School 
for Advanced Studies (SISSA), Via Beirut 2-4, 1-34014 Trieste, Italy 
(Dated: February 11, 2005) 

Using exact diagonalizations and Green's function Monte Carlo simulations, we have studied 
the zero-temperature properties of the quantum dimer model on the triangular lattice on clusters 
with up to 588 sites. A detailed comparison of the properties in different topological sectors as 
a function of the cluster size and for different cluster shapes has allowed us to identify different 
phases, to show explicitly the presence of topological degeneracy in a phase close to the Rokhsar- 
Kivelson point, and to understand finite-size effects inside this phase. The nature of the various 
phases has been further investigated by calculating dimer-dimer correlation functions. The present 
results confirm and complement the phase diagram proposed by Moessner and Sondhi on the basis 
of finite-temperature simulations [Phys. Rev. Lett. 86, 1881 (2001)] . 



I. INTRODUCTION 

The investigation of spin-liquid phases is currently a 
very active field of research, partly - but not only - be- 
cause of their possible connection to the superconductiv- 
ity observed in several cuprates. The definition of a "spin 
liquid" is itself a matter of debate. Following the work 
of Shastry and Sutherland on a two-dimensional model 
whose exact ground state is a product of dimer singletSj 1 
the word is sometimes used to designate phases in which 
the spin-spin correlation function decays exponentially 
fast with distance at zero temperature. However, such 
phases often exhibit other types of long-range order, like 
dimer order, which manifest themselves as non-decaying 
correlation functions involving more than two spins^ In 
that respect, the word liquid is not appropriate, and it 
should arguably be reserved to systems in which all cor- 
relation functions decay exponentially fast at large dis- 
tance. This discussion would be quite academic if the 
only characteristic of such liquids was the absence of any 
kind of order, but following the pioneering work of Wen^ 
it is well admitted by now that such liquids can exhibit 
another property known as topological order. In the ther- 
modynamic limit, the ground state (when defined on a 
topologically nontrivial domain) exhibits a degeneracy 
not related to any symmetry and referred to as topo- 
logical degeneracy. These degenerate ground states live 
in topological sectors which cannot be connected by any 
local operator. 

The realization of such phases in quantum spin mod- 
els is still preliminary though. The best candidates are 
frustrated magnets for which quantum fluctuations are 
known to destroy magnetic long-range order, but their 
ground-state properties are very difficult to access, and 
when definite conclusions are reached, it is usually be- 
cause the presence of some kind of long-range order 
(dimer, plaquette, etc.) can be established** The main 
difficulty is in a sense technical: A good diagnosis would 
require to study large enough clusters, but this is not 
possible since quantum Monte Carlo simulations of frus- 



trated antiferromagnets are plagued with a very severe 
minus sign problem. 

In that respect, effective models such as the quantum 
dimer model (QDM) are extremely interesting. Although 
their relationship to actual Heisenberg antiferromagnets 
is not a simple issued they describe resonance processes 
typical of strongly fluctuating frustrated quantum mag- 
nets while retaining the possibility to be analyzed by 
standard techniques such as quantum Monte Carlo. This 
possibility was first exploited by Moessner and Sondhi, 6 
who developed a finite-temperature Monte Carlo algo- 
rithm to study the QDM on a triangular lattice defined 
by the Hamiltonian: 

H = vJ2(\~)(~\ + \//)(//\) 

where the sum runs over all plaquettes including the 
three possible orientations. The kinetic term controlled 
by the amplitude t changes the dimer covering of ev- 
ery flippable plaquette, i.e., of every plaquette contain- 
ing two dimers facing each other, while the potential 
term controlled by the interaction v describes a repul- 
sion (v > 0) or an attraction (v < 0) between dimers 
facing each other. Since a positive v favors configura- 
tions without flippable plaquettes while a negative v fa- 
vors configurations with the largest possible number of 
flippable plaquettes, the so-called maximally flippable 
plaquette configurations (MFPC), one might expect a 
phase transition between two phases as a function oiv/t. 
The actual situation is far richer though. As shown by 
Moessner and Sondhi, who calculated the temperature 
dependence of the structure factor, there are four differ- 
ent phases (see Fig.^l: i) A staggered phase for v/t > 1, 
in which the ground-state manifold consists of all non- 
flippable configurations; ii) A columnar ordered phase 
for v/t sufficiently negative; iii) An ordered phase adja- 
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FIG. 1: Phase diagram of the QDM on the triangular lattice 
as a function of v/t after Ref. lu- 



cent to it called the \/12 x V12 phase where the struc- 
ture factor develops a low-temperature peak at an in- 
termediate value of the wave-vector and which probably 
consists of resonating plaquettes which make a 12-site 
unit-cell pattern; iv) A liquid phase with a featureless 
and temperature-independent structure factor. This last 
phase has been interpreted as a short-range resonating 
valence bond (RVB) phase in which all correlations de- 
cay exponentially, even at low temperature^ It is sepa- 
rated from the staggered phase by a special point, the 
Rokhsar-Kivelson (RK) point (v/t = 1). At this point, 
the ground-state manifold contains all staggered configu- 
rations plus another configuration of zero energy, the sum 
of all possible configurations. In that particular state, 
the dimer-dimer correlation function has been shown to 
decay exponentially with distance at zero temperature^ 

While the arguments put forward by Moessner and 
Sondhi in favor of this phase diagram are quite convinc- 
ing, this proposal calls for further investigation for several 
reasons. First of all, this zero-temperature phase dia- 
gram was inferred from finite-temperature results, and a 
direct analysis of ground-state properties would be wel- 
come. Besides, the location of the phase boundaries is to 
a large extent unknown, as well as the nature of the quan- 
tum phase transitions between the different phases. Fi- 
nally, and more importantly, a direct investigation of the 
RVB phase would help to understand its physical proper- 
ties. In particular, this phase is expected to be liquid in 
the strong sense of the word, and, as such, to have topo- 
logical degeneracy, but this could not be checked by the 
finite temperature calculation of Moessner and Sondhi. 
In fact, this property, which lies at the root of Ioffe et al's 
proposal for q-bits^ has never been directly observed. As 
noticed by Ioffe and collaborators on the basis of exact 
diagonalizations of small clusters, the finite-size effects 
are still huge for the small cluster sizes accessible with 
that technique, and no conclusion regarding the thermo- 
dynamic limit could be reached. 

In this paper, we address all these issues, and answer 
most of them, by a careful investigation of the zero- 
temperature properties of the model essentially based on 
the implementation of a Green's function Monte Carlo 
(GFMC) algorithm. The huge finite-size effects are 
shown to be a natural consequence of the subtle inter- 
play between the cluster geometry, the order parameter 
of the underlying phase (if any) and the topological sec- 



tor. This analysis turns the finite-size effects into a very 
powerful tool to investigate the ground-state properties 
of the model as a function of v/t. In particular, we have 
obtained strong evidence in favor of topological degen- 
eracy close to the RK point, and we have been able to 
locate the transition between the two ordered phases with 
a reasonable accuracy. 

The paper is organized as follows. In section II, we 
present the minimal technical background necessary to 
understand the results presented in the following sec- 
tions. To keep this technical section as small as possi- 
ble, some details have been relegated into Appendices. 
We then discuss the various phases and the transitions 
between them in Section III, and finally we conclude in 
Section IV. 



II. TECHNICAL BACKGROUND 

A. Topological sectors 

On the triangular lattice, the QDM has conserved 
quantities defining topological sectors (TS)&2ii£ On fi- 
nite clusters, these conserved quantities are defined by 
the parity of the number of dimers intersecting a given 
line which satisfies two properties: 1) It is closed, or it 
ends at the boundaries of the cluster; 2) It does not divide 
the cluster into two disconnected pieces. The first condi- 
tion ensures that the parity is conserved when applying 
the Hamiltonian, while the second condition guarantees 
the possibility to construct configurations with both par- 
ities. On a cylinder, the only choice is a line going from 
one end to the other, while on a torus, one can choose 
any closed loop that goes around one of the axis of the 
torus. Since the sectors defined by two lines that can be 
deformed into each other continuously are the same, one 
ends up with two sectors for a cylinder and four sectors 
for a torus. In the following, we will work exclusively 
with clusters defined on a torus. 



B. Finite clusters 

On the triangular lattice, it is possible to construct 
two types of clusters which keep all the symmetries of 
the infinite latticeiii In terms of the basis vectors Ui and 
u 2 [with ui = (1,0) and u 2 = (1/2,73/2)], they are 
defined by two vectors: 



Ti = lui + mu 2 

T 2 = -mui + (I + ra)u 2 



(1) 

(2) 



with I or m = for type-A clusters (see Fig. |5J) and 
I = m for type-B clusters (see Fig. |3| . The clusters have 
the geometry of tori obtained by identifying sites of the 
infinite lattice modulo the vectors Ti and T 2 . 

The number of sites N of such clusters is given by the 
simple formula: N = I 2 + Irn + m 2 , leading to N = I 2 
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FIG. 2: Example of a type-A cluster with 16 sites. 
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FIG. 4: Degeneracy of the TS as a function of the parity of 
1/2. Three out of four TS are degenerate since they contain 
the same configurations rotated by an angle n/3. 
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TABLE I: Topological sectors in which the various states 
and phases discussed in the text can be found. The symbol * 
means that the size of the cluster must be a multiple of 12 to 
be compatible with the \/T2 x \/l2 phase. 



FIG. 3: Example of type-B clusters with 12 sites showing the 
two different choices of unit cell. 



sites for type-A clusters and N = 3 x l 2 sites for type- 
B clusters. Since the number of sites must be even to 
accommodate dimers, we consider only clusters with I 
even. 



C. Topological sectors, order, and cluster size 

The four topological sectors that we denote with (0,0), 
(0,1), (1,0), and (1,1) can be defined by considering loops 
along Ui and 112, where the first (resp. second) number 
refers to the loop along Ui (resp. u 2 ), and (respec- 
tively 1) stands for even (respectively odd). Remarkably, 
even on a finite cluster, three of the four sectors are al- 
ways equivalent: This property is a consequence of the 
correspondence between sectors under rotations (see Ap- 
pendix[SJ and applies equally to type-A and type-B clus- 
ters, and, therefore, the three sectors have the same spec- 
trum. However, which sectors are degenerate depends on 
the size of the cluster, according to the rule of Fig. 0] 

If the system is in a completely liquid phase with 
no long-range order at all, all topological sectors are 
expected to be equivalent up to finite-size effects, and 
each topological sector should contain one replica of the 
ground state. However, if some dimer order develops, 
this is no longer true, and the ground state will only be 
found in some sectors. The rule depends on the cluster 



size and on the type of order. The columnar order can 
only appear in the (0,0) sector, and, according to the rule 
of Fig. ^] in the (0,1) and (1,0) sector when 1/2 is odd. 
On the other hand, the y/12 x y/l2 phase is not frus- 
trated in clusters with a number of sites multiple of 12, 
i.e., all type-B clusters and type-A clusters with I mul- 
tiple of 3. Further, it is only compatible with the (0,0) 
sector for 1/2 even and with the (1,1) sector when 1/2 is 
odd. This rule is again equally applicable to both type-A 
and type-B clusters. The above properties are summa- 
rized in Table [I] They will prove helpful to identify the 
phases and their boundaries since we only have access to 
ground-state properties in different topological sectors. 



D. Numerical methods 

Since the size of the Hilbert space grows exponentially, 
exact diagonalizations could only be performed for one 
cluster of type B (12 sites) and two clusters of type A 
(16 and 36 sites). A few typical results are listed in the 
Appendix [5] However, on the basis of these results only, 
not much could be said regarding the issues of topological 
degeneracy in the thermodynamic limit and of the phase 
transitions. 

Fortunately, since all non-diagonal matrix elements are 
non-positive, the zero-temperature GFMC method could 
be implemented H In particular, we make use of the algo- 
rithm with a fixed number of walkers described in Rcf. El 
The algorithm is ergodic at sufficiently large v/t, except 
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for the staggered configurations, which anyway play no 
role in the parameter range v/t < 1 in which we are in- 
terested. However, the large energy barriers present for 
negative values of v/t, typically v/t < —1, make it dif- 
ficult to browse the configuration space efficiently. In 
particular, for very negative v/t, the probability to sam- 
ple states around the MFPC starting from a generic one 
becomes exponentially small, introducing serious ergod- 
icity problems in the algorithm (Appendix 0. 

Results have been obtained with this algorithm for 
clusters with up to 588 sites. While no guiding func- 
tion was required very close to the RK point, a guiding 
function favouring clusters according to the number of 
flippable plaquettes had to be implemented to keep a 
good statistics when decreasing v/t. 



III. THE PHASE DIAGRAM 




FIG. 5: Topological gap for type-A clusters as a function of 
the linear size I = N 1 ^ 2 of the cluster. The error bars are 
smaller than the symbol size. 



We now turn to the presentation of the results we 
have obtained starting from the RK point and decreas- 
ing toward the most negative values of v /t affordable by 
GFMC. As mentioned previously, the analysis is based 
on a comparison of the ground-state energy in differ- 
ent topological sectors. Since the sectors (0,1) and (1,0) 
are always degenerate with (0,0) or (1,1), depending on 
the parity of 1/2, a comparison of these last two sectors 
is sufficient. In the following, we will denote by 2?(o,o) 
[resp. i)] the total ground-state energy in sector (0,0) 
[resp. (1,1)] and by AE the difference between them: 
AE = 2?(i,i) — 2?(o,o), which we call the topological gap. 
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A. The RVB phase 

In this section we present our numerical results close 
to the RK point, where a disordered spin- liquid phase is 
expected. We, therefore, look for evidences of topological 
degeneracy in the thermodynamic limit by a size-scaling 
study. Our analysis has shown that finite-size effects are 
very different for clusters of type A and B, and we prefer 
to discuss them separately. 



1. Type-A clusters 

We have calculated the topological gap for type-A clus- 
ters with up to the 22 x 22 cluster and the results for the 
topological gap are shown in Fig. [5] The finite-size effects 
are very strong, and using only the exact diagonalization 
data, which extend up to I = 6, it would be hopeless 
to draw any conclusion about the thermodynamic limit. 
Instead, with the sizes reachable within GFMC, a clear 
tendency appears: The topological gap seems to decrease 
with the cluster size with an oscillatory behavior. It has 
been shown in Ref. ^3 that the same behavior must be 
observed in the correlation function of auxiliary fermions 



FIG. 6: Oscillatory fits of the topological gap per site near 
the RK point for type-A clusters. The plot is in a log-linear 
scale. 

introduced in the Pfaffian method at the RK point. This 
correlation function behaves along Ui (or u.2) as&iS 

G(R) cx R- 1/2 exp(-i?/£) , R -> 00, (3) 

where £ has both a real and an imaginary part given 
by ~ 0.83 + i0.127Tiii Accordingly, we have tried to 
fit our results for the topological gap with an oscillatory 
function defined as: 

Ae(Z) = ^ e -H-[r l ]i cos (lm[r ^ + ) , (4 ) 

where Ae = AE/N is the energy difference per site. 
Fig. [fj] displays these fits for three different ratios of v/t 
near the RK point. As expected, the fit is very accurate 
for v/t close to the RK point, but becomes worse as v/t 
decreases. 

The real and imaginary parts of the inverse correlation 
length extracted from these fits are reported in Fig.[7| 
The continuous lines correspond to the values of v/t for 
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FIG. 7: Real and imaginary part of the inverse correlation 
length £ as a function of v/t. At the RK point, the values 
quoted are the analytical result of Ref. HH For type-A clus- 
ters, both real and imaginary parts of are shown. The 
continuous lines connect the point for which the fit has a satis- 
factory accuracy (the error bars are of the order of the symbol 
size) . 



which the fit with Eq. Q was suficiently reliable. As 
v/t — > 1, these values connect smoothly to the analytic 
results obtained by Ioselevich and collaborators at the 
RK point. 

Because Re[£ _1 ] remains finite away from the RK 
point, we can conclude that the topological gap vanishes 
exponentially in the thermodynamic limit in a finite pa- 
rameter range close to the RK point, as expected in the 
RVB phase. For v/t < 0.85, Im[£ _1 ] saturates at the 
value corresponding to the yT2 x vT2 crystal phase and 
we shall comment on this fact in the section UlI Bl The 
non-trivial behavior of the topological gap for different 
sizes of the cluster suggests to look at clusters of differ- 
ent geometry. Therefore, we continue our analysis with 
type-B clusters. 



2. Type-B clusters 

Now, we consider the type-B clusters of size N = 3 x I 2 
with I up to 14. The finite-size effects of the topolog- 
ical gap for this kind of cluster are much simpler than 
for type-A clusters. Indeed, AE changes sign accord- 
ing to (— l)'/ 2 , and its absolute value smoothly decreases 
with the cluster size without any oscillations. Follow- 
ing Ref. 0, we may relate the topological gap in type- 
B clusters in the close vicinity of the RK point to the 
fcrmionic Green's function in the direction Ui + 112 ( "B 
direction" in the terminology of Ref. \ldS) . Because of 
the special symmetry properties along this direction, the 
corresponding correlation length is purely real at the RK 
point and its value is ~ 0.76. An accurate calcula- 
tion of the fermionic Green's function at the RK point 
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FIG. 8: Fits of the topological gap per site near the RK point 
for type-B clusters. We choose to multiply the total energy 
by the prefactor N~ 7 ^ 8 to fit with a pure exponential function 
[see Eq. ©]. 



along the lines of Ref. gives 

G(R) cx R- 1/4 exp(-R/£) , R^-oo. (5) 

Thus in the vicinity of the RK point we may fit the topo- 
logical gap with^ 
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const(-l)'/ 2 V3 i/£ 

Jl/4 



(6) 



where the correlation length £ is real (v3 I is the length of 
the vectors Ti and T2 in type-B clusters). Indeed, even 
away from the RK point, the "B direction" conserves its 
special symmetry and we expect to see no incommen- 
surate oscillations in Ae(l). Our numerical calculations 
confirm this behavior, and as for type-A clusters, the fits 
are good close to v/t = 1 if one only takes into account 
the sizes much larger than the correlation length (see 
Fig-EJ- This exponential decay confirms the presence of 
topological degeneracy close to the RK point. 

Remarkably the (— \) 1 / 2 oscillations in Ae(l) connect 
smoothly the RVB phase to the ^fY2 x Vl2 crystal 
phase, where the ground state must alternate between 
the (0,0) and (1,1) topological sectors. This indicates 
that short-range VT2 x VT2 correlations are probably al- 
ready present in the ground state of the RVB phase. 



3. Correlation function 

To demonstrate the absence of long-range correlations 
in the vicinity of the RK point, we have also calculated 
the dimer-dimer correlation function (D(x)D(0)}, where 
D(x) is the dimer operator which is equal to if there 
is no dimer on bond x and to 1 if there is one. In the 
RVB liquid phase, we expect all dimer orientations to 
be equally probable far enough from the reference dimer 
-D(O). Since the average dimer concentration is 1/6, 
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FIG. 9: Correlation function of two parallel dimers as a func- 
tion of their separation x along ui for v/t = 0.97. Inset: 
Dimer-dimer correlation function on the 108-site cluster for 
v/t = 0.97 in the topological sector (1,1). Dashed (resp. solid) 
links correspond to probabilities smaller (resp. larger) than 
1/36, while the thickness of the links is proportional to the 
absolute value of the distance to 1/36. The reference dimer 
is in the lower-left corner. 






FIG. 10: A simplified scheme of the the V12 x V12 phase. 
The unit cell (shaded hexagon) consists of 12 sites and con- 
tains three resonating rhombi plaquettes. This construction 
neglects correlations and resonances between different plaque- 
ttes of the same unit cell and between different unit cells. 



(D(x)D(0)) must tend to 1/36 as x — > oo. This liquid 
behavior of the dimer-dimer correlations has been con- 
firmed at the RK point by Moessner and Sondhi, 6 «& who 
have proven that, at v/t = 1, (D(x)D(Q)) tends to 1/36 
exponentially with increasing distance x. 

Not surprisingly, we find the same behaviour for v/t 
close to 1, which lends further support to the existence 
of an RVB phase away from the RK point. An example 
of short-range dimer correlations at v/t = 0.97 is shown 
in Fig. El 



B. The transition between the RVB and the 
intermediate (\/l2 x \/V2) phases 

Upon decreasing v/t, we expect a phase transition from 
the short-range RVB phase to an ordered phase with 
dimer or plaquette order. According to Moessner and 
Sondhi, this should occur somewhere between v/t = 
and v/t = 1, presumably around v/t = 0.7, and the phase 
on the other side of the transition should be the ordered 



7 12 x vT2 phase (see Fig. ITU1 for a schematic picture of 
the order or Refs. 1(1 1 til for a more detailed discussion) . 

From the results obtained on type-A clusters, one may 
suspect that the phase transition takes place around 
v/t = 0.85. Indeed, around this value Im[£ -1 ] saturates 
at a value close to ir/6, and at the same time the pro- 
posed fit Q becomes inaccurate. Note that ir/6 is ex- 
actly the value one would expect for the periodicity of 
the topological gap in type-A clusters. Indeed, the pe- 
riod of the \/l2 x yl2 crystal in the directions of Ui 
and U2 is 6 lattice spacings, and, therefore, the period- 
icity of the topological gap is 12 lattice spacings, since 
the ground-state topological sector changes as / changes 
by 6. The saturation value of Im[^ _1 ] reported in Fig.[7| 
is somewhat larger than 7r/6, which is related to the in- 
applicability of the fitting formula J3J for small clusters. 
As one can measure directly from the last half-period of 
oscillations in Fig. EI at v/t < 0.85 the correct value of 
the oscillation wave vector Im[£ -1 ] in fact equals tt/6 to a 
very good precision, in perfect agreement with yfl2 x \/T2 
ordering. 

Another estimate of the phase transition point comes 
from the analysis of type-B clusters. The number of sites 
in these clusters is always a multiple of 12, therefore, they 
can accomodate a vT2 x VT2 crystal without defects. As 
we have already mentioned, the sign oscillations of the 
topological gap AE = (— l)'/ 2 |Ai?| agree both with the 
RVB liquid and with the VT2 x \/12 crystal phase, and 
it is only the size dependence of \AE\ that may indicate 
the phase transition. 

We have plotted in Fig. El the absolute value of the 
topological gap as a function of 1/yfN for several val- 
ues of v/t ranging from 1 to negative values. There is a 
clear change of behavior for large clusters. For v/t close 
to 1, the topological gap monotonously decreases with 
the size, as discussed earlier, while for v/t close to 0, the 
topological gap seems to diverge with the size. This is 
actually what we expect in the \/l2 x vT2 phase. Indeed, 
since this order can only appear in one topological sec- 
tor (see section llTC")) . considering the wrong topological 
sector will require a one-dimensional defect, with an en- 
ergy cost proportional to y/N. We could not reach large 
enough cluster to actually check this scaling form, but 
the data are consistent with it. The precise location of 
the transition is also difficult to pin down, but the change 
of curvature for the large sizes occurs around v/t = 0.7, 
which may be taken as a candidate value for the transi- 
tion point. 

To summarize, from the analysis of topological gaps in 
type-A and type-B clusters, the phase transition point 
between the RVB and the VT2 x \/l2 phases may be es- 
timated between v/t ~ 0.7 and v/t ~ 0.85. However, 
our results do not shed much light on the nature of the 
phase transition (see, e.g., a possible transition scenario 
with 0(4) symmetry in Ref. Il6|) . We may only remark 
that the correlation lengths reported in Fig. [7| admit the 
possibility of a second-order phase transition: the cor- 
relation lengths for type-A and type-B clusters merge 
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FIG. 11: Absolute value of the topological gap for type-B 
clusters as a function of l/\/JV for several values oiv/t. There 
is a clear change of behavior between v/t — 1 and v/t = 0. 
Lines are guides for the eye. 



together quickly as v jt decreases from the RK point (the 
correlation length should become isotropic and divergent 
at the second-order transition point). Such a transition 
would also explain our difficulties in locating the exact 
transition value of v/t from simulations on small clusters. 
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FIG. 12: Upper panel: ground state dimer-dimer correlations 
of the 108-site cluster for v/t — in the topological sector 
(1,1). The same plotting scheme as in the inset of Fig. [5] 
Lower panel: the schematic plot of dimer-dimer correlations in 
the classical state of Fig. 1101 Thick lines represent peaks of the 
correlation function higher than 50% above the average value 
1/36. The gray dimer in the bottom left corner represents the 
reference link. 



C. The intermediate (VT2 x VT2) phase 

In addition to the analysis of the topological gaps in 
type-A and type-B clusters, we may obtain a more direct 
evidence of the existence of the VT2 x %/l2 phase from the 
dimer-dimer correlation functions. Upon decreasing v/t 
from the RK point, the long-distance correlations grad- 
ually increase. The resulting pattern for 108 sites and 
v/t = is shown in Fig. 1121 

We can identify the structure of these correlations as 
belonging to the intermediate \/l2 x \JY2 phase. For- 
tius purpose, we simulate the dimer-dimer correlations 
in a classical state corresponding to the schematic repre- 
sentation of the intermediate phase in Fig. ^| Namely, 
we assume that all dimer configurations within a fixed 
Vl2 x vT2 crystal of hexagons come with equal prob- 
ability. The resulting pattern of increased dimer-dimer 
correlations is presented in the lower panel of Fig. ^] 
Even though this construction is oversimplified, we ex- 
pect that it correctly predicts the main features of the 
correlation function. A direct comparison with the up- 
per panel of Fig. ^] indicates a y/l2 x \/l2 correlation 
pattern at v jt = 0. 



D. The transition between the intermediate and 
the columnar phases 

According to Moessner and Sondhi, a phase transition 
to the columnar phase is expected when v /t becomes suf- 
ficiently negative, although they were not able to locate 
it precisely. If we look at appropriate clusters, this transi- 
tion should correspond to a change of ground-state topo- 
logical sector. Indeed, while in the columnar phase the 
ground state always appears in the (0,0) sector [and, for 
1/2 odd also in (0,1) and (1,0)], in the \/l2 x \/l2 phase, 
it appears in the (0,0) sector if Z/2 is even and in the (1,1) 
if 1/2 is odd (see section Hi C|) . Therefore, if we look at a 
cluster multiple of 12 with 1/2 odd, the transition from 
the vT2 x \/T2 to the columnar phase should correspond 
to a change of ground-state topological sector from (1,1) 
to (0,0). In other words, the topological gap is expected 
to vanish at that point. To avoid mixing the results of 
clusters with different geometries, we have concentrated 
on available type-B clusters with appropriate geometries, 
namely clusters with 12, 108, 300 and 588 sites. The re- 
sults are reported on Fig. El F° r 12 sites, they are ob- 
tained by exact diagonalizations, and the level crossing 
occurs at v/t = —2.022. For the other sizes, they have 
been obtained with GFMC. Because of ergodicity prob- 
lems, appearing at negative v/t (see Appendix , we 
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FIG. 13: Absolute value of the topological gap as a function 
of v/t for type-B clusters with 12, 108, 300 and 588 sites. The 
dashed lines are extrapolations to locate the level crossing 
between the (0,0) and (1,1) topological sectors. Inset: Linear 
extrapolation of the level crossing points. 



could not reach the level crossing, but we could get close 
enough to it to allow a meaningful extrapolation of the 
results. The results obtained are reported as a function 
of the inverse linear size 1 / y/N in Fig. ^] It seems rea- 
sonable to approximate this size dependence by a linear 
function, which leads to an estimate v/t = —0.75 ± 0.25 
for the transition point between the \/l2 x vL2 phase 
and the columnar phase. 



E. The columnar phase 

Finally, although there is little doubt concerning the 
nature of the phase for large and negative v/t, we have 
tried to calculate the dimer-dimer correlation function in 
this phase as well for completeness. As mentioned before, 
there is an ergodicity problem, and we cannot reach that 
phase with GFMC. However, for clusters for which this 
transition appears as a level crossing between topologi- 
cal clusters, we can expect to get an idea of the correla- 
tions taking place in this phase by calculating them for 
a value of v/t as close as possible to the transition and 
in the topological sector compatible with the columnar 
phase. We have thus calculated the dimer-dimer corre- 
lation function for the 108-site cluster for v/t = —0.8 
in the topological sector (0,0). The results are plotted in 
Fig. 1141 A clear columnar pattern appears, in contrast to 
the correlations in the intermediate phase (see Fig. I12[). 
Note that one columnar configuration has been selected 
from the others. This is due to the fact that the rota- 
tional symmetry is broken when one considers a single 
topological sector. 




FIG. 14: Dimer-dimer correlation for the 108-site cluster at 
v/t = —0.8 in the topological sector (0,0). A clear columnar 
pattern appears. The same plotting scheme as in the inset of 
Fig. El 



IV. CONCLUSIONS 

In this paper, we have performed a thorough analysis 
of the ground-state properties of the QDM on the tri- 
angular lattice with exact diagonalizations and GFMC. 
New results with respect to previous investigations could 
be obtained thanks to the combination of two factors: i) 
The large sizes reachable with GFMC (up to 588 sites); 
ii) The precise identification of the interplay between the 
cluster geometry, the topological sector and the nature 
of the ground state. 

The most important result of this paper is to provide 
the first direct numerical evidence that the four topolog- 
ical sectors on a torus are indeed degenerate in the ther- 
modynamic limit in a finite range below the RK point, 
thus confirming the identification of this phase as an RVB 
phase by Moessner and Sondhi. The finite-size effects to- 
ward this unique ground-state energy as revealed by the 
behaviour of large clusters agree, in the vicinity of the 
RK point, with the analytic perturbative expressions of 
Ioselevich and collaborators>i& Furthermore, fitting with 
their perturbative formula allowed a meaningful extrac- 
tion of the correlation length also away from v/t = 1. 

The other phases of the phase diagram of Moessner 
and Sondhi, namely the \/l2 x \f\2 and the columnar 
phases, have also been confirmed by an investigation of 
the dimer-dimer correlation function, and a careful anal- 
ysis of the crossing between two topological sectors for 
certain clusters has allowed us to conclude that the phase 
transition takes place at v/t = —0.75 ± 0.25. 

The transition between the RVB phase and the Vl2 x 
\/l2 phase is more difficult to pin down, but the behavior 
of the topological gap for large clusters suggests that it 
takes place between v/t = 0.7 and v/t = 0.85. 

The main issue to be resolved is the nature of the phase 
transition between the RVB phase and the Vl2 x vT2 
phase. The correlation length increases away from the 
RK point, but we could not obtain a reliable evidence 
that it diverges, as would be the case for a continuous 




TS(0,0) 2 TS(0,0) 



FIG. 15: Three cut-lines labelled by 1, 2 and 3 on a type A 
cluster with (a) 16 sites (1/2 even) and (b) 36 sites (1/2 odd). 



phase transition. In that respect, it would be very useful 
to implement a dynamical version of the GFMC algo- 
rithm to calculate the excitation gaps to vison and non- 
vison excitations inside a given topological sector (thus 
extending beyond the RK point the results of Ref . ITtI) . 
Work is in progress along these lines. 
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APPENDIX A: RELATIONS BETWEEN 
TOPOLOGICAL SECTORS 

As stated in section III CI and summarized in Fig. 0] 
the topological sectors (0,1) and (1,0) are always degen- 
erate either with (0,0) or (1,1) depending on the cluster 
size. The reason for this is simply that three out of four 
topological sectors are related by 7r/3 rotations. In or- 
der to know to which topological sector (0,1) and (1,0) 
are connected to, it is useful to keep track of the parity 
of the number of dimers intersecting a line in the third 
direction, i.e., 112 — Ui (see Fig. 11511 . A rotation by ir/3 
corresponds to a cyclic permutation of the parities along 
the three cut lines. For the sector (0,0), the parity along 
the third line is clearly if 1/2 is even and 1 if Z /2 is odd, 
as illustrated in Fig. ^Jfor type- A clusters. Thus, if 1/2 
is even, the (0,0) sector is left unchanged by a rotation, 
and (0,1) and (1,0) are connected to (1,1). Instead, if 1/2 
is odd, the (0,0) sector is connected to (0,1) and (1,0), 
and (1,1) is left unchanged. 
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FIG. 16: Examples of maximally flippable plaquette config- 
urations in all sectors for 1/2 even. 



APPENDIX B: EXACT DIAGONALIZATIONS 

For further reference, we list in Table [H] the energies 
obtained with exact diagonalizations for the 4x4 and 
the 6x6 clusters of type A and for the 12-site cluster of 
type B. 



APPENDIX C: ERGODICITY AND MAXIMALLY 
FLIPPABLE PLAQUETTE CONFIGURATIONS 

If t = and v < 0, the ground-state manifold is 
highly degenerate: It contains all the maximally flippable 
plaquette configurations. These configurations have the 
largest possible number of flippable plaquettes, i.e., N/2. 
For 1/2 odd, they can only be found in the sectors (0,0), 
(0,1) and (1,0). For 1/2 even, such configurations are 
present in all sectors (see Fig. ^| for an example) . For 
large negative but finite v/t, the degeneracy is lifted at 
second order in t/v in favour of columnar states^ The 
investigation of the ground-state properties by GFMC 
in this range of v/t suffers from the complicated energy 
landscape of the phase space. The large energy bar- 
riers that separate the configurations tend to trap the 
GFMC simulations in metastable minima, and exponen- 
tially long simulation times would be necessary to con- 
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v/t 


TS 


12 sites 


16 sites 


36 sites 


0.8 


00 


-0.69505 


-1.02863 


-2.10806 




11 


-0.84755 


-0.96885 


-2.20081 


0.6 


00 


-1.48237 


-2.10582 


-4.36730 




11 


-1.70960 


-2.02804 


-4.55215 


0.4 


00 


-2.31590 


-3.21999 


-6.74713 




11 


-2.58514 


-3.13239 


-6.98904 


0.2 


00 


-3.17522 


-4.36297 


-9.21881 




11 


-3.47308 


-4.26664 


-9.48840 


0.0 


00 


-4.05317 


-5.52971 


-11.76017 




11 


-4.37228 


-5.42488 


-12.03779 


-0.2 


00 


-4.94745 


-6.71735 


-14.35767 




11 


-5.28161 


-6.60465 


-14.62961 


-0.4 


00 


-5.85765 


-7.92454 


-17.00373 




11 


-6.20000 


-7.80513 


-17.25894 


-0.6 


00 


-6.78415 


-9.15091 


-19.69534 




11 


-7.12645 


-9.02644 


-19.92244 


-0.8 


00 


-7.72760 


-10.39677 


-22.43356 




11 


-8.06006 


-10.26927 


-22.61782 



verge to the real minimum. This situation is different 
from the usual behavior where the convergency to the 
ground state is related to the physical gap. To make sure 
that such problems were absent in the results presented 
in this paper, we have systematically compared the out- 
come of several simulations, starting with different pools 
of walkers. For each cluster size, this has defined the 
lowest reachable value oi v/t. 



TABLE II: Total ground-state energy obtained with exact 
diagonalizations for the topological sectors (0,0) and (1,1) and 
for different sizes. 
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